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Cn I Abstract. I review our current understanding of the properties of strongly 

H ' interacting matter at high temperatures, based upon numerical calculations in lattice 

QCD. I discuss the chiral and deconfining aspects of the QCD transition, the equation 
of state, fluctuations of conserved charges, color screening, meson correlation functions, 

CO ■ and the determination of some transport coefficients. 
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1. Introduction 

r-| ■ It is expected that strongly interacting matter shows qualitatively new behavior at 

temperatures and/or densities which are comparable or larger than the typical hadronic 
scale. It has been argued by Hagedorn that production of hadronic resonances may 
O ■ Ificid to a limiting temperature above which hadronic matter can no longer exist [1]. 

^ . This limiting temperature was re-interpreted by Cabibbo and Parisi as the transition 

VO ■ temperature to a new state of matter that consist of quarks [2j. Indeed, based of 

cn . asymptotic freedom one would expect that relevant degrees of freedom at very high 

^ . temperatures and/or densities are quarks and gluons which are no longer subject to 

confinement [3j, i.e. deconfined. Since the quarks and gluons interact weakly at high 
temperatures the deconfined matter is analogous to the plasmas and is also called the 
/S . quark gluon plasma (QGP) [^. In particular, it is characterized by color screening that 

c^ . is similar to the well- know Debye screening |5]. 

The existence of such deconfining transition was first shown using strong coupling 
expansion in lattice QCD P E], followed by numerical Monte-Carlo studies of the lattice 
SU{2) gauge theory which confirmed it [8l[9],[T0]. The deconfinement of gluons was seen 
as onset of color screening [HI [9] and rapid increase in the energy density that indicate 
liberation of many new degrees of freedom [10] . The onset of color screening and rapid 
increase in number of degrees of freedom remain the key signatures of deconfinement 
even today. 

Since these pioneering studies QCD at finite temperature became quite a large 
sub-field of lattice QCD. One of the obvious reasons for this is that the transition from 
hadronic to partonic matter cannot be described in perturbation theory. But even at 
high temperatures the physics is non-perturbative beyond the length scales l/{g^ T) 
[g"^ being the gauge coupling) pAj. Therefore lattice QCD remains the only ab-initio 
tool for theoretical understanding of the properties of strongly interacting matter under 
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extreme condition which is important for the physics of the early universe [121 [131 [S] 
and heavy ion colhsions |15j . 

Symmetries play an important role in understanding the transition to QGP. For 
zero quark masses the QCD Lagrangian has an approximate chiral symmetry S'f/^(3) 
as well as the axial symmetry Ua{^)- The later is broken by quantum effects, while 
the former is spontaneously broken by the QCD vacuum |T6]. At high temperatures 
the chiral symmetry is expected to be restored |T7]. Therefore there should be a chiral 
transition in QCD. The nature of the chiral transition was first discussed in Ref. [T8] . 
where it was argued that for three light quark flavors the transition should be first order, 
while for two flavors it could be second order belonging to the 0(4) universality class. 
Furthermore, the axial symmetry is also effectively restored at high temperatures as 
the breaking of the axial symmetry is controlled by instanton density [H] that vanishes 
at infinite temperature [17j. If the breaking of axial symmetry at the chiral transition 
temperature is reduced the chiral transition in the two flavor theory could be first order 
[18] . I will comment on this issue in section |5l 

In the opposite limit of infinitely heavy quark masses the QCD partition function 
has Z{3) symmetry [9]. This symmetry is broken at high temperatures due to color 
screening. In this case we have a deconfining transition which is first order. The 
Z{3) symmetry is a good symmetry for sufficiently large quark mass. For intermediate 
quark masses the transition is a crossover. The regions of the first order transitions 
are separated from the crossover region by lines of second order transition that belong 
the Z{2) universality class. The region of the first order deconfinement transition was 
mapped out in Refs. [HI |20l [21]. The endpoints of this first order region correspond to 
a quark masses of the order of the charm quark mass, niq ~ 1.4GeV [HI [20] . 

The boundary of the first order chiral transition was studied using effective models, 
see Refs. [22] [23] for recent works (also references therein). It was found that for 
three degenerate quark fiavors the endpoint of the region where the transition is first 
order corresponds to a pseudo-scalar meson mass of (110 ± 40)MeV. There is no direct 
evidence for the first order chiral transition for three very light quark fiavors so far in 
lattice QCD. This maybe due to lattice artifacts. The existing calculations that use 
staggered fermion formulation give only an upper bound for the pseudo-scalar meson 
mass smaller than 70MeV [211 [23 [2E], which is not inconsistent with effective model 
analysis mentioned above. The end-point of the first order transition, or in more general 
case the line that separates the first order transition from the crossover region belongs 
to the Z{2) universality class. This also means that for some value of the strange quark 
mass, rn^^^ we have a tricritical point that connects Z{2) and 0(4) lines as well as the 
first order region. Calculations in linear sigma model suggest that m^^^ corresponds 
to the kaon mass of about (1700 — 1850)MeV [23]. On the other hand, as we will see, 
lattice QCD calculations indicate a much smaller value of rn^'~'^ . The above picture 
of the finite temperature QCD transition is summarized in Fig. [1] (also known as the 
Columbia plot). 

For the physical quark masses the transition is likely to be an analytic crossover 
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Figure 1. The QCD transition as function of the quark masses. 



(for earlier works see Refs. P5|I2U] ). Therefore, the discussion of the deconfinement and 
the chiral aspects of the QCD transition at finite temperature and the determination 
of the corresponding temperature requires a special care. In fact, it is only the chiral 
aspects of the transition that allow a determination of the corresponding transition 
temperature. In sections |3] and H] I discuss the deconfinement aspects of the QCD 
transition in terms of chromo-electric screening and fluctuations of conserved charges. 
The chiral aspects of the QCD transition and the determination of the corresponding 
transition temperature are then discussed in section [5l Equation of state is one of the 
most important characteristic of hot strongly interacting medium. It will be discussed in 
sectionlH Further insight into the properties of medium can be gained by studying meson 
correlation functions and spectral functions. They are also important for heavy ion 
phenomenology. Therefore section [7] is dedicated to the discussion of meson correlation 
functions and spectral functions. In the next section I will review some basic information 
on lattice QCD calculations which is relevant for the present discussions. 



2. Basics of lattice QCD 



To study QCD non-perturbatively we use lattice gauge theory [30]. In this formalism 
a field theory is defined in a gauge-invariant way on a discrete space-time domain. 
This serves at least two purposes: a) to provide an ultra-violet cut-off for the theory, 
restricting highest momentum to vr/a (a being the lattice spacing), and b) to evaluate the 
path integrals in the Euclidean formulation stochastically using importance sampling. 
On the lattice the fundamental degrees of freedom of a theory with local SU{?)) 
gauge symmetry are fermion fields tpx that reside on the sites of the lattice and carry 
flavor, color and Dirac indeces, which we suppress through the most of this paper, and 
bosonic, gauge degrees of freedom that in the form of S'f/(3) matrices f/^_^ reside on 
links. Sites on a four- dimensional lattice are labeled with x = (r, x). The link variables 
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are related to the usual gauge fields: U^^^ = exp{igaAf^{x)), where a is the lattice spacing 
and g is the gauge coupling. 

The theory is defined by the partition function 



Z = DUDiljDijexp{-S) (1) 

where the action 

S = S, + Sf (2) 

contains gauge, Sg and fermionic, Sf parts. The latter part is bi-linear in fields and has 
the form 

Sf = J2^,D,ijg (3) 

where Dg = D + rrig is the fermion matrix and the sum goes over the quark fiavors. In 
the simplest formulation the lattice gauge action can be written as 

^5 = /3 E (l - l^^ Tr U^,,U,^f.,uUU^Ul) , (4) 

(3 = Q/g"^ with g"^ being the bare gauge coupling. This is the Wilson gauge action [30] . 
It is easy to see that expanding the link variables in a the above expression gives the 
well known Yang-Mills action X^x o!^^ / ^{FavY ^P ^o corrections that are suppressed as 
a^ . Thus calculations with Wilson gauge action have discretization errors of order a^. 
These discretization errors can be removed by the Symanzik improvement program that 
relies on higher order difference scheme, i.e. considering combination of link variables 
that extend beyond the elementary square [3T| . 

Th expectation value of an operator O is then given by 

(O) = ^ f DUD^DtpOexpi-S). (5) 

Integration over the fermion fields (which are Grasmann variables) can be carried out 
explicitly: 

Z = jDU{detDg[U]rexp{-Sg) ^ J DU exp{-S,ff), (6) 

where Se/f = Sg — Uf In det Dq[U] is the effective action. The fermion determinant 
det Dq[U] describes the vacuum polarization effects due to the dynamical quarks and 
makes the effective action non-local in gauge variables. For this reason simulations 
with dynamical quarks are very resource demanding and the quenched approximation 
is often employed, where det Dq[U] is set to 1. This approximation will be relevant for 
the discussion of the meson spectral functions. 

The most commonly used fermion formulations include : the Wilson fermion 
formulation, staggered fermion formulation [32] and the domain wall fermion (DWF) 
formulation [331 EH ES] . The Wilson formulation breaks the chiral symmetry of the QCD 
Lagrangian. As the result quark masses acquire an additive renormalization which makes 
the calculations quite complicated and numerically expensive. The domain wall fermions 
preserve all the symmetry of the QCD Lagrangian at a price on introducing the extra 
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fifth dimension. Tlierefore tlie computational costs associated witli DWF formulation 
are proportional to the extent of the fifth dimension and are typically 30 to 100 times 
larger than the ones with Wilson fermions. The staggered formulation describes four 
fermion flavors (also called tastes) in the continuum limit. The full Ul{4) x Uji{4) 
symmetry of the four flavor theory is reduced to Us{l) x Ub{1), where Us{l) is the 
subgroup of the SUa{4:) symmetry. Thus the staggered formulation preserves a part of 
the chiral symmetry. This makes finite temperature calculations with staggered fermions 
attractive for two reasons. First, due to the remnant chiral symmetry there is no additive 
renormalization of the quark mass. This makes staggered fermions computationally 
inexpensive. Second, the remnant chiral symmetry makes it easy to study the chiral 
aspects of the finite temperature QCD transitions. To study arbitrary number of flavors 
with the staggered fermion formulation the rooting trick is used. This amounts to 
replacing Uf in Eq. ([6]) with n//4. The validity of this procedure was studied in details 
in recent years [361 EH [Ml EH SQl SH US] though some open issues remain (see e.g. Ref. 
[3H1I13] ). 

The standard staggered fermion formulations has discretization errors of order a^. 
They come either from the lattice effects in the fermion dispersion relation and formally 
start at order g^a'^, and from the breaking of taste symmetry, which starts at order 
g^a"^. The former lead to large cutoff effects in the thermodynamic quantities at high 
temperatures. The later lead to the distortion of the hadron spectrum at non-zero lattice 
spacing. In particular, out of 16 pseudo-scalar meson masses only one vanishes in the 
zero quark mass limit. The remaining 15 pseudo-scalar meson have masses proportional 
to a ( the squared mass of these mesons goes like g^a"^ for sufficiently small lattice 
spacings) and thus vanish only in the continuum limit. These discretization errors can 
be eliminated by Symanzik improvement program. Using 3-link terms in the lattice 
Dirac operator one can eliminate tree level a^ errors and largely reduce cutoff effects in 
thermodynamics quantities |11] . Such improvement of the staggered fermion formulation 
is mandatory to control discretization effects at high temperature [H]. To reduce errors 
due to taste symmetry breaking so-called fat links are used. These are combinations of 
the usual link U^^^ and products of links along different paths. One of the possibilities 
is to construct fat links by using APE smearing [45j : a link variable Ux,fj. is replaced by 
a weighted average of itself and a sum of the 3-link paths connecting the same sites as 

U.,, ^ f/;^ = (1 - Qc)Ux,, + c E U.,,Ux+o,,Ul^f.,.. (7) 

Here c is a weight factor, < c < 1. More complicated paths can be used as well. 
In particular, using combination of different paths up to length seven it is possible to 
eliminate taste breaking effects up to order g^a"^ |16]. The fat links described above 
are not elements of the gauge group SU{3). While this is not a problem in principle, 
it tuns out that projecting the fat links to elements of SU{3) or f/(3) gauge group 
further reduces the breaking of taste symmetry [171 HH |19]. The staggered fermion 
actions that used in large numerical calculations take certain combination of fat links 
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and improvement of the quark dispersion relation and go under acronym pA, asqtad and 
HISQ/tree. In HISQ/tree action the fat hnks are projected to U{3) [50]. The stout 
action used by Budapest- Wuppertal collaboration only includes fat links [5ll |52l [271 EH] . 
To evaluate the path integral ([5]) stochastically, an ensemble of Nu gauge 
configurations, weighted with exp(— Sg//), is generated using Monte Carlo methods. 
The state of the art Monte-Carlo algorithm that can handle arbitrary number of quark 
flavors and most often used in QCD thermodynamics goes under the name of the rational 
hybrid Monte-Carlo (RHMC) algorithm [5l]- The expectation value of the operator is 
then approximated by the ensemble average: 

1 Nu 

(o)^TrEo.(f/), (8) 

^ 1=1 

where Oi{U) is the value of the operator O, calculated on i-th configuration. This 
method works as long as Seff is real. When quark chemical potentials are included, 
e.g. as in calculations at non-zero net baryon density, 5*/ and Se/f become complex and 
Monte-Carlo methods do not work. This is the well known sign problem. 

All the quantities calculated in lattice simulations are dimensionless, i.e. are 
calculated in units of the lattice spacing a. The lattice spacing is varied by varying 
the lattice gauge coupling /3 ~ 1/(7^. To set the scale in lattice simulations we need to 
calculate a given quantity and assume that it has the known physical value. Different 
choices are possible and the lattice spacing determined by different quantities should 
be the same up to discretization errors that vanish when the continuum limit is taken. 
Common choices used in the literature are the scale parameters tq and vi determined 
by the static quark anti-quark potential V{r): 

r^ — !,=,„ = 1.65, r^ — 1,=,^ = 1.00 (9) 

Since the static quark anti-quark potential cannot be measured experimentally the 
value of the parameters r^ and r\ in physical units in the continuum should be 
determined using some other experimentally measured quantity. The most precise 
determination of r\ comes from using the pion decay constant /^ as an input, which gives 
r\ = 0.3106(20)fm [55] • Furthermore, studying the shape of the static potential from 
the above value of ri one can obtain the value of tq = 0.468(4)fm [56]. Alternatively one 
can use fx to set the lattice spacing. While this procedure is more straightforward, for 
the staggered fermion formulation it is also more problematic as one expects that due 
to taste symmetry breaking fx has larger discretization errors than the static potential. 
The above discussion referred to the zero temperature case. To consider the case 
of non-zero temperature, T the Euclidean time extent has to be 1/T and periodic and 
anti-periodic boundary condition have to be imposed on the boson and fermion fields 
[57] . Thus the temperature is related to the lattice spacing and the temporal extent 
of the lattice, T = l/(iV^a). Taking the continuum limit at fixed temperature implies 
Nr ^ CO while keeping the aspect ratio N^/Nr fixed, where A^'^- is the spatial extent of 
the lattice. 
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3. Deconfinement : color screening 

3.1. Color electric screening and order parameter for deconfinement 

As discussed in the introduction the deconfining transition in SU{N) gauge theories 
is a true phase transition related to Z(N) symmetry |||. The order parameters of this 
phase transition are the expectation value of the Polyakov loop and the Polyakov loop 
correlator 

L(T) = {-TrW{x)), W{x) = H Uo{t,x), (10) 

^^ r=0 

CpL{r,T) = i^(TrP^(r)Trl^(0)). (11) 

The Polyakov loop transforms non-trivially under Z{N) transformation and the 
expectation value of the Polyakov loop is zero in the confined phase and non-zero above 
the transition temperature. The correlation function of the Polyakov loop is related 
to the free energy of a static quark anti-quark pair [9]. The qualitative change in the 
behavior of the Polyakov loop and its correlator above the phase transition temperature 
is related to color screening. As it was pointed out already in Ref. |9] the Polyakov 
loop and the Polyakov loop correlator require renormalization to be interpreted as the 
free energy of an isolated static quark or the free energy of a static quark anti-quark 
iQQ) pair. More precisely, they are related to the free energy difference of a system 
with static QQ pair at some temperature and the system without QQ pair at the same 
temperature. Since in the zero temperature limit the free energy of a static quark anti- 
quark pair should coincide with the static potential the Polyakov loop renormalization 
is determined by the normalization constant of the static potential, namely 

UUT) = exp(-c/(2r))L(T) = exp(-Foo(T)/(2r)) (12) 

CpL{r,T)= exp{-F{r,T)/T + c/T), F^{T) = \imF{r,T), (13) 

where c is an additive normalization that ensures that the static potential has a certain 
value at a given distance. In the above expressions we made explicit that the free 
energy of an isolated static quark anti-quark pair is half the free energy of QQ pair 
at infinite separation. In the confined phase the free energy of a static quark anti- 
quark pair is proportional to cr(T)r at large distances r, as expected. The effective 
temperature dependent string tension cr(T) is non-zero below the phase transition 
temperature [611 [62] . Consequently the free energy of an isolated static quark is infinity 
and LreniT) =0. In the deconfined phase Foo(T) is finite due to color screening. In 
particular at high temperatures F^o = — 4asm£)/3 at leading order, where rriD = gT 
is the leading order Debye mass. Here and in what follows we consider only the case 
A^ = 3. The next-to-leading order correction for F^o has been calculated in Refs. [63|[6^ 

I The deconfinement transition for other gauge groups have been discussed in Refs.j SHilS^ISU] . 
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and was found to be small. For the free energy of a static QQ pair at distance r we 
have l65l 



la' 
F(r,T) = --^exp(-2mz5r). (14) 

This result is in contrast to the free energy of static charges in QED which is 

—aexp{—mDr)/r, and is the consequence of non-Abelian nature of interactions. A 

static quark anti-quark pair could be either in a singlet or in an octet state, and thus 

in a fixed gauge) one can define the so-called color singlet and octet free energy [651 Efi] 



exp(-Fi(r,T)/T + c/T) = i(Tr[Lt(x)L(y)]), (15) 

exp{-Fs{r,T)/T + c/T) = l(TrLt(a;)TrL(2/)) 

o 

-^(Tr[Lt(x)L(i/)]). (16) 

Then the Polyakov loop correlator can be written as the thermal average over the singlet 
and the octet contributions P |65l [66] 

Cpi(r,T) = ^exp(-Fi(r,T)/T + c/T) + ^exp(-F8(r,T)/T + c/T).(17) 

The singlet and octet free energies can be calculated at high temperature in leading 
order HTL approximation [71], resulting in 

F^{r,T) = - -^expi-mnr) ^, (18) 

6 r 6 

Fa{r,T) = -—exp{-mDr) — . (19) 

or 6 

The singlet and octet free energies are gauge independent at this order. Inserting the 

above expressions into Eq. fITTj) and expanding in a^ one can see that the leading order 

contributions from the singlet and octet channels cancel each other and one recovers 

the leading order result for F{r,T) in Eq. (TT71) . It was not clear how to generalize the 

decomposition of the free energy into the singlet and octet contributions beyond leading 

order. Recently using the effective theory approach, namely the potential non-relativistic 

QCD (pNRQCD) at finite temperature [721 [73] it was shown that this decomposition 

indeed holds at short distances [64J- As we will see below the singlet free energy turns 

out to be useful when studying screening numerically in the lattice calculations. 

The QCD partition function does not have the Z{3) symmetry. This symmetry 

is broken by the quark contribution. In physical terms this means that a static quark 

anti-quark pair can be screened already in the vacuum by light dynamical quarks. The 

free energy does not rise linearly with distance but saturates at some distance, i.e. 

§ The terms color singlet free energy and octet free energy are misleading as only F{r,T) has the 
meaning of the free energy of static quark anti-quark pair. Furthermore, at low temperatures both 
Fi and _Fg are determined by color singlet asymptotic states |67j . The color screening of static quark 
quark interaction was discussed in Ref. [68j , while color screening of charges in higher representations 
was studied in Refs. [691I70I. 
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Figure 2. The Polyakov loop as function of the temperature in 2+1 flavor QCD 
[74l [56] and in pure gauge theory [62l [75] . 



we see string breaking. For light dynamical quarks the corresponding free energy F^o 
is twice the binding energy of static-light meson and thus is proportional to Aqcd- 
Therefore, the Polyakov loop is expected to be 0{1) in the transition region. At very 
high temperatures, however, the screening in the QCD and purely gluonic plasma is 
expected to be similar. The above discussion of the free energy of static quark anti- 
quark pair still holds and the only difference is in the value of the Debye mass, which 
now becomes rriD = gTJl + Nf/6. Here Nj being the number of light quark flavors, 
i.e. Nf = 3. 

In Fig. |2] I show the recent lattice results for the renormalized Polyakov loop 
in 2+1 flavor QCD with physical quark masses [71] compared to the results in pure 
gluodynamics [621 CS]- The temperature scale in pure gauge theory was set using the 
value y/a = 465.2MeV for the zero temperature string tension. We see a relative smooth 
increase in the Polyakov loop in 2+1 flavor QCD starting at temperatures of about 
160MeV. Earlier results for the renormalized Polyakov loop in 2 and 3 flavor QCD with 
signiflcantly heavier quark masses have been reported in Refs. [TBI [77| The behavior of 
the Polyakov loop in pure glue theory and full QCD is quite different in the transition 
region, while at high temperatures it is qualitatively similar. 

The free energy of a static QQ pair as well as the singlet free energy was calculated 
using p4 action on 16^ x 4 and 24^ x 6 lattices for almost physical light quark 
masses, namely mi = rris/lO [78l [79]. There are also calculations with Wilson fermion 
formulation with much heavier quark masses [80] as well as preliminary results for the 
stout action [81]. Qualitatively these results are similar to the ones obtained with pA 
action. A very detailed calculation of the free energy of QQ pair in pure glue theory 
was presented in Refs. [611 |62l [82] , while the singlet free energy in Coulomb gauge was 
studied in Refs. [751 183]. 

In Fig. |3]I show the numerical results for the free energy and the singlet free energy 
obtained with p4 action on 16^ x 4 lattice. At low temperatures the qualitative behavior 
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of the free energy and singlet free energy is similar. The singlet free energy agrees with 
the zero temperature potential shown in the figure as the black line, while the free energy 
differs from it by the trivial color factor Tin 9. At high temperatures the free energy 
and the singlet free energy show very different behavior, in particular the temperature 
dependence is much stronger in F{r,T). The singlet free energy agrees with the zero 
temperature potential at short distances, while this is not the case for the free energy 
even when the factor Tin 9 is taken into account. This behavior of the QQ free energy 
at high temperatures is consistent with the partial cancellation of the color singlet and 
color octet contributions as described above. Similar difference between the free energy 
and singlet free energy have been observed in the pure glue theory |6Tl [75l IH2] • 

Color screening means that r(Ti(r, T)—Foo{T)) as well as r{F{r, T)— Too(T)) should 
decay exponentially at large distances. Therefore in Fig. H] I show these combinations. 
At distance rT > 0.8 we indeed see the expected exponential falloff. Fitting the lattice 
data for r{Fi{r, T) — F^{T)) to an exponential one can determine the Debye mass non- 
perturbatively. This was done for pure glue theory and 2+1 flavor QCD. It turns out 
that the temperature dependence of the Debye mass is well described by the leading 
order result but its value is about factor 1.4 larger both in pure glue theory and 2+1 
flavor QCD [79j [jjj. This means that the dependence of the Debye mass on quark flavors 
is well described by perturbation theory. The Debye mass is non-perturbative beyond 
leading order [HU [HH] . The non-perturbative enhancement of the Debye mass over the 
leading order result is expected based on the gap equation calculation of the Debye mass 
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Figure 3. The free energy (left) and the singlet free energy (right) of a static QQ pair 
calculated in 2+1 flavor QCD on 16^ x 4 lattices at different temperatures [75]. The 
solid black line is the parametrization of the zero temperature potential calculated in 
Ref. [87]. 



The qualitative features of the singlet free energy discussed above are not specific to 
the Coulomb gauge. In fact, the singlet free energy defined in terms of Wilson loops show 
very similar behavior to the one calculated in the Coulomb gauge [88]. Furthermore, 



II The Debye mass obtained in the calculations with Wilson fermions is somewhat larger [80] 
calculations are needed to resolve this discrepancy. 
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the Debye mass can also be defined from tlie long distance beliavior of tlie electric gluon 
propagator [89l |90l EH |92l [931 [M] • Calculations in SU{2) and SU{3) gauge theories show 
no gauge dependence of the extracted Debye mass within statistical errors [921 193] . The 
extracted screening masses are in agreement with the ones obtained from the singlet 
free energy 



3.2. Color magnetic screening and dimensional reduction 

Contrary to ordinary plasmas that have no magnetic screening the static chromo- 
magnetic fields are screened in QGP. This is due to the fact that unlike photons gluons 
interact with each other (the stress tensor is non-linear in QCD). Magnetic screening is 
non-perturbative, i.e. it does not appear at any finite order of pertubation theory and is 
needed to cure the infrared divergences of the pressure rendering it finite but also non- 
perturbative at order g^ [TT]. In analogy with electric screening magnetic screening can 
be studied using gauge fixed magnetic gluon propagators. However, the large distance 
behavior of the magnetic propagators is more complex and much more susceptible to 
finite volume effects. In Landau type gauges as well as in Coulomb gauge the magnetic 
propagator shows oscillating behavior instead of decaying exponentially [95]. Thus no 
magnetic mass can be defined. However, one can define a common exponentially falling 
envelope for different gauges that then gives a screening mass of about 0.5g^T [95] . 
A magnetic mass of similar size was found in analytic approach of Refs. [961 EZl EH] . 
Thus the corresponding effective magnetic screening length is larger than the electric 
screening length as one would expect for weakly coupled QGP l5ll. 

Spatial Wilson loops can also be used to define a length scale associated with 
static magnetic fields. They obey area law at any temperature governed by spatial 
string tension as{T). The numerical results for the spatial string tension calculated 
in 2+1 flavor QCD with p4 action and mi = rris/lO are shown in Fig. [5l At low 



% Speculations on existence of a magnetic screening length of similar size have been presented in Ref. 
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Figure 5. (Left)The string tension in units of rp as function of the temperature 
calculated on A',- = 4, 6 and 8 lattices. For better visualization we scale the 
temperature axis by To = 200MeV. (Right) T/^/ojT) calculated on A^r = 4, 6 and 
8 lattices compared with the prediction of dimensional reduction indicated by the 
line. The uncertainty in the prediction of the dimensionally reduced theory is show by 
dashed lines. 



temperatures crs(T) coincides with the zero temperature string tension, while at high 



temperatures Jas{T) ~ T. This behavior can be understood in terms of dimensionally 
reduced effective theory. At high temperatures there is a separation of different scales 
2nT ^ rriD ^ g^T. If we integrate out the modes associated with the scale 27rT the 
thermodynamics and the screening lengths can be described in terms of an effective 
field theory, which at the lowest order is the 3 dimensional adjoint Higgs model in the 
confined phase ^Qi [TOU [102] 



^3" 



t/^x— Tr F,,(x)Fi,(x) 

93 



Tr[A,Ao(x)]2 
+ m^TrAo(x)2 + AA(TrAo(x)2)2 . (20) 

The static electric fields become adjoint scalars in the three dimensional effective theory. 
The couplings of the effective theory can be calculated in perturbation theory at any 
order |10H 1102] . The spatial string tension should be identified with the string tension 
of this effective three dimensional theory a^. The gauge coupling constant in three 
dimensions g^ has the dimension of the mass, and at the lowest order (tree level ) 
g^ = g^T. On dimensional grounds a^^ = cg^ with c being a constant. This explains the 
linear rise of the string tension at high temperatures. The coefficient c was calculated 
in the three dimensional adjoint Higgs model |103[ 1104] . The three dimensional gauge 
gs coupling was calculated to 2- loop |105] . Combining this with the known value of c 
we get the weak coupling prediction for the spatial string tension which is compared 
with the lattice data in Fig. [51 Quite surprisingly, the weak coupling result for the 
spatial string tension works well down temperatures of 200MeV, where the effective 
field theory is not supposed to work. Note that in the conformal strongly coupled gauge 
theory T/J(7s{T) ~ 1/A^/^ (see e.g. Ref. |106] and references therein). Thus the 
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ratio T/J(Js{T) is temperature independent and vanishes in the hmit of infinite t'Hooft 
couphng, A — )■ oo. This is in sharp contrast with lattice calculations. To deal with 
this problem one can modify the AdS metric to break the conformal invariance and 
get a QCD-like theory. This approach is known as AdS/QCD. Attempts to model the 
temperature dependence of the spatial string tension using AdS/QCD were discussed 
inRef. [TnTllTnSllTnQ] . 

3.3. The spectrum of gauge invariant screening masses 

There are several screening masses that can be studied in the framework of the 3 
dimensional effective theory |110[ IllH I112[ I103[ 1104] as well as in the original four 
dimensional theory |113[ I114[ I115[ [62] . Comparison of the screening masses calculated in 
these two theories can provide information on the validity of the dimensionally reduced 
effective theory. The screening masses can be classified according to J^^ quantum 
numbers. Here R is the analog of C-parity. It is related to the symmetry with respect 
to the following transformation of the static electric field, Aq — > —Aq. All the operators 
transform trivially under charge conjugation in the adjoint Higgs model |103] . In the 0+"'" 
channel there are two distinct screening masses: screening mass m{Al) corresponding 
to the "electric" correlator of Tr^Q, and the screening mass mo corresponding to the 
"magnetic" correlator TrF^^-. These correlators can mix but this mixing is suppressed 
for weak coupling. The lightest screening mass tue in the channel was proposed 
as the non-perturbative definition of the Debye mass |116] since the negative R parity 
prevents mixing with the magnetic sector. For asymptotically small couplings we expect 
rriG < rriE < m{AQ) and the large distance behavior of the Polyakov loop correlators is 
actually governed by the "magnetic" screening mass mo |116| 1117] . For the interesting 
temperature range, however, m{Al) turns out to be the smallest and m^ turns out to be 
the largest of the above screening masses |103[ 1104] . Furtheremore, there is almost no 
mixing between the "electric" and "magnetic" correlators in the O"*"^ channel |103l 1104] . 
Thus the long distance behavior of the Polyakov loop is determined by the screening 
mass corresponding to TtAq correlator ( see also discussion in Ref. |117] ). In Fig. [6] I 
show the screening masses extracted from the Polyakov loop correlators in SU{2) and 
SU{3) gauge theories and compared to the corresponding screening masses calculated in 
the effective three dimensional theory. There is fairly good agreement with the screening 
masses calculated in the full theory and in the effective theory. 

The lightest screening mass in the O"*"^ channel was calculated in perturbation 
theory at next-to-leading order |118] and the corresponding result is also shown in 
the figure. The next-to-leading order result significantly overestimates the screening 
masses, while the leading order result seems to work fine. The next-to-leading order 
calculations also explains to some extent why the screening mass in the TtAq channel 
is the smallest jllSj . An interesting observation has been made in Ref. |115] . where 
screening masses corresponding to the correlators of the real and imaginary parts of the 
Polyakov loop have been calculated. In terms of the effective 3d theory these correlators 
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Figure 6. The lightest screening mass in 0++ channel calculated from the Polyakov 
loop correlators in SU{2) gauge theory [62l I114i 1119) (left) and SU{?>) gauge theory 
[6n I115J (right). Also shown as open symbols are the results of the calculations in 
the 3d effective theory jl03| I104J . The solid and dashed lines are the perturbative 
predictions at leading order and next-to- leading order [118) . 



correspond to the correlators of TrAg and TrAjj respectively. The ratio of the screening 
masses extracted from these correlators should be 2/3 at leading order and the lattice 
calculations of Ref. |115] confirm this result. 

4. Deconfinement : Fluctuations of conserved charges 



As mentioned before due to the infamous sign problem lattice QCD Monte-Carlo 
simulations are not possible at non-zero quark chemical potentials. However, the 
pressure at non-zero chemical potentials can be evaluated using Taylor expansion. The 
Taylor expansion can be set up in terms of the quark chemical potentials /x^, /i^ and 
/is, or in terms of the chemical potentials corresponding to baryon number B, electric 
charge Q and strangeness S of hadrons 



p_ 
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^{^^B|Ty^{^^Q/Ty^{^^s/TY^ ^^^^ 

While Taylor expansion can be used to study the physics at non-zero baryon density 
the expansion coefficients are interesting on their own right as they are related to the 
fluctuations and correlations of conserved charges. The diagonal expansion coefficients 
are related to second and higher order fluctuations of conserved charges 



X2 



VT^ 
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xf = :^((iV^)-3(iVi)^) (25) 

X^ = ^, {{N'x) - 15(iV^)(iVi) + 30(iVi)3) , (26) 

while off-diagonal expansion coefficients are related to correlations among conserved 
charges, e.g. 

Xfi'' = p^(iVxiV^). (27) 

The quark chemical potentials are related to the chemical potential of baryon number, 
electric charge and strangeness 
1 2 

1 1 

11 

f^S = g/^B - g/iQ - f^S- (28) 

Therefore the expansion coefficients in quark chemical potential Xuds ^^^ related to the 
hadronic ones Xbqs- I^ particular, for the second order expansion coefficients we have 

X2 =1 (X2 +Xi + X2 + 2Xu + 2Xn + 2x^1 

'' ■ ^' '' " ■ us I n^ ,ds A^,ud 



X2 = g [^Xl + Xt + X2- 4Xii + 2x?? - 4xi? 

X2 ~ X25 



Xi? = \ {2X^2 - X2 - X2 + Xii - 2xil + X?0 

xfi^ = -\{x'2 + xTi + xil) , 



X?i = ^ (X2 - 2X^1^ + Xil) . (29) 

Fluctuations and correlations of conserved charges are sensitive probes of 
deconfinement. This is because fluctuation of conserved charges are sensitive to the 
underlying degrees of freedom which could be hadronic or partonic. Fluctuations of 
conserved charges have been studied using different staggered actions [281 ESI [ZH I120[ 
fT2n fT22l fT23l [TM fT25l fT26l [T271 [T28] . As an example in Fig. 1 show the quadratic 
strangeness fluctuation as the function of the temperature calculated with HISQ /tree 
and stout actions. Fluctuations are suppressed at low temperatures because conserved 
charges are carried by massive strange hadrons (mostly kaons). They are well described 
by Hadron Resonance Gas (HRG) model at low temperatures. Strangeness fluctuations 
rapidly grow in the transition region of T = (160 — 200) MeV as consequence of 
deconfinement. At temperatures T > 250 MeV strangeness fluctuations are close to 
unity that corresponds to the massless ideal quark gas. Similar picture can seen in 
other fluctuations, in particular for light quark number fluctuation, xu which is also 
shown in Fig. [71 The more rapid rise of xi i^ the transition region is a quark mass 
effect. In fact, similar quark mass dependence is seen in the HRG model. At high 
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temperatures the strange quark mass plays little role. In fact the difference between 
the strange and the light quark number susceptibilities is zero within errors. There is a 
good agreement between the calculations performed with stout and HISQ/tree action 
as the continuum limit is approached. Let me finally note that strangeness fluctuations 
have been calculated very recently with Wilson action for the physical values of the 
quark masses, and these calculations confirm the staggered result |129] . 
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Figure 7. Strange quark number (left) and light quark number fluctuations (right) 
calculated with stout [127] and HISQ/tree actions |i56.. The lattice results are 
compared with the prediction of the HRG model shown as a black line. 
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Figure 8. Quark number correlation x^*! (left) and Cbs (right) calculated with stout 
[127] and HISQ/tree actions |1261 1128) . The lattice results are compared with the 
prediction of the HRG model shown as a black line. 

Correlations of conserved charges, in particular quark number correlations are 
also sensitive to the relevant degrees of freedom. At infinitely high temperatures the 
correlations between different quark numbers should vanish due to the weakly interacting 
nature of QGP. In the high temperature region quark number correlations have been 
calculated in perturbation theory |130[ I13H 1132] and in the framework of dimensionally 
reduced effective theory |133j . The deviation from the ideal gas limit turned out to be 
quite small. At low temperatures, where hadrons are the relevant degrees of freedom 
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such correlations naturally arise as different quark ffavors are bound in hadrons. In Fig. 
[8] I show the results of the lattice calculations performed with stout and HI SQ /tree 
action for Xil- These results demonstrate the expected features: at high temperatures 
the correlations are very small while at low temperatures they are compatible with 
HRG results. The fact that Xn has a minimum is due to interplay between meson and 
baryon contributions. Mesons and baryons contribute to Xn with opposite signs. At low 
temperatures the negative contribution from mesons dominates Xil- As the temperature 
increases baryon contribution becomes significant and eventually the dominant one 
because of the larger density of states in the excited baryon sector. It is also interesting 
to consider the following combination 

Cbs = -3xf//x2 = 1 + ^i^^' (30) 

that was first introduced in Ref. |134j . The Boltzmann suppression at low temperature 
is canceled out in this combination and at high temperatures it should approach one. 
The lattice results obtained with stout and HISQ/tree actions are also shown in Fig. ID 
Again, at low temperatures the lattice results are well described by HRG model, while 
at high temperatures Cbs is close to one as expected for weakly interacting quarks. 

Higher order expansion coefficients have been studied in detail for p4 action and 
light quark mass mi = O.lnis [121] . In Fig. [9] I show the fourth and sixth order 
fluctuations of the electric charge. The forth order coefficient develops a peak in the 
transition region and at high temperatures approaches the ideal quark gas value. The 
sixth order expansion coefficient has a peak around the transition temperature and 
becomes negative above the transition temperature. It approaches zero from below as 
the temperature increases as expected in the weakly interacting QGP. The six order 
expansion coefficient vanishes in the ideal gas limit as well as in the next order of 
perturbation theory [130j . Based on the discussion discussion of the second order 
expansion coefficients we would expect that at low temperatures also the higher order 
expansion coefficients should be reasonably well described by HRG. However, this is not 
the case. As one can see from Fig. [9] the lattice results fall below the HRG prediction. 
This is due to the large cutoff effects related to taste symmetry breaking in the p4 
action. Taste symmetry breaking effects lead to distortions of the hadron spectrum. 
If these distortions are taken into account in the HRG calculations a good agreement 
with the lattice data can be achieved |135l 11361 1137] . The HRG model predictions with 
the distorted spectrum are shown in Fig. [9] as the dashed and dotted lines and are in 
reasonably good agreement with the lattice data. Preliminary results for the higher 
order expansion coefficients have been recently obtained with HISQ/tree action [124j . 
The qualitative features of these expansion coefficients remain the same, however, the 
location of the peak is shifted to smaller temperatures. The observed shift in the peak 
position is consistent with the shift in the chiral transition temperature observed in 
the calculations with HISQ/tree action [56]. The higher order expansion coefficients 
are sensitive to the singular part of the free energy density in the vicinity of the chiral 
transition temperature. Therefore I will come back to the discussion of these quantities 
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Figure 9. The fourth (left) and sixth (right) Taylor expansion coefficients for the 
electric charge chemical potential calculated for pA action [121] . The prediction of the 
HRG model are shown as solid black line. The dotted and dahsed lines correspond to 
the HRG model with the distorted hadron spectrum |136[ 1137) . 



in section |5l which is dedicated to the chiral aspects of the QCD transition. 

In the high temperature region the fourth order expansion coefficient was also 
calculated with pA action on Nr = 4:, 6 and 8 lattices |123] . While no continuum 
extrapolation was performed the numerical data suggest that at high temperatures the 
fourth order expansion coefficient will approach the ideal quark gas value from below 
|123] . In the same study charm quark number fluctuations have been considered using 
partially quenched approximation, i.e. no charm quark loops. It was found that the 
ideal quark gas is a good approximation also in this case |123] . 

Finally let us compared the lattice results for quark number susceptibilities at high 
temperatures with the prediction of resummed perturation theory |130| 1131] . In Fig. 
[To] I show the lattice results for light quark number susceptibilities obtained with p4, 
stout and HISQ/tree actions compared with the resummed perturbative results. Since 
the difference between light and strange quark number susceptibilities is negligible for 
T > 400MeV in the figure I also show the strange quark number susceptibilities obtained 
with asqtad action [28]. The continuum extrapolated stout results agree well with HTL 
resummed perturbative result jl30]. Whereas, the A^,- = 8 asqtad results, that are 
expected to be close to the continuum limit, agree with the results obtained using the 
next-to-leading log (NLA) approximation |131] . Clearly more lattice calculations in the 
high temperature limit are requited to settle this issue. It is also interesting to compare 
the results of the lattice calculations with the results obtained in strongly coupled gauge 
theories using AdS/CFT correspondence. The result from AdS/CFT calculations shown 
in the figure as the solid black line is significantly below the lattice results cJ. 

+ The conserved charges considered in these calculations are not exactly the quark numbers, see 
discussion in Ref. [138]. 
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5. Chiral transition 

The Lagrangian of QCD has an approximate S't/^lS) chiral symmetry. This symmetry 
is spontaneously broken in the vacuum. The chiral symmetry breaking is signaled by 
non-zero expectation value of the quark or chiral condensate, (ijjip) 7^ in the massless 
limit. This symmetry is expected to be restored at high temperatures and the quark 
condensate vanishes. There is an explicit breaking of the chiral symmetry by the small 
value of u, d and s quark masses. While due to the relatively large strange quark 
mass {ms — 100 MeV) S'f/4(3) may not be a very good symmetry its subgroup S'[/a(2) 
remains a very good symmetry and is relevant for the discussion of the finite temperature 
transition in QCD. If the relevant symmetry is 5'[/a(2) the chiral transition is expected 
to be second order for massless light {u and d) quarks belonging to the 0(4) universality 
class [18]. Recent calculations with pA action confirm this pictures |139] . In other words 



m 



phys 



> m 



TCP 



contrary to calculations in effective linear sigma model. This also 
means that for non-zero light quark masses the transition must be a crossover. The 
later fact seems to be supported by calculations in Ref. [27]. The Ua{^) symmetry 
is explicitly broken in the vacuum by the anomaly but it is expected to be effectively 
restored at high temperatures as the non-perturbative vacuum fluctuations responsible 
for its breaking are suppressed at high temperatures. If the Ua{^) symmetry is restored 
at the same temperature as the SUa{'2) symmetry the transition could be first order 
|18] . Recent calculations with staggered as well as with domain wall fermions suggest 
that Ua{^) symmetry gets restored at temperature that is significantly higher than the 
chiral transition temperature |140[ 1141] . 

For massless quarks the chiral condensate vanishes at the critical temperature 
T° and it is the order parameter. Therefore in the lattice studies one calculates the 
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chiral condensate and its derivative with respect to the quark mass, called the chiral 
susceptibility. For the staggered fermion formulation most commonly used in the lattice 
calculations these quantities can be written as follows: 

(V^^).,^ = iiv^Tr(D,-i), (31) 

Xm,q{T) = Uf ^ = Xq4isc + Xq,con q = l,S, (32) 

where the subscript x = t and x = denote the expectation value at finite and zero 
temperature, respectively. Furthermore, Dg = niq ■ 1 + D is the fermion matrix in the 
canonical normalization and nj = 2 and 1 for light and strange quarks, respectively. 
In Eq. (13 2 p we made explicit that chiral susceptibility is the sum of connected and 
disconnected Feynman diagrams. The disconnected and connected contributions can be 
written as 

X.,...c = Y^^ {{{TrD;^) - {TrD^r} , (33) 

Xq,can = - ^TiJ^i D^^x, 0)D;\0, x) ) , q = l,s. (34) 

The disconnected part of the light quark susceptibility describes the fluctuations in the 
light quark condensate and is analogous to the fluctuations of the order parameter of 
an 0{N) spin model. The second term {Xq,con) arises from the explicit quark mass 
dependence of the chiral condensate and is the expectation value of the volume integral 
of the correlation function of the taste non-singlet scalar operator ipip. Let me note 
that in the massless limit only xi,disc diverges. In the next subsections I will discuss the 
temperature dependence of the chiral condensate and the chiral susceptibility as well as 
the role of universal scaling in the transition region. 

5.1. The temperature dependence of the chiral condensate 

The chiral condensate needs a multiplicative, and also an additive renormalization if 
the quark mass is non-zero. Therefore the subtracted chiral condensate is considered 

In Fig. [11] I show the results for {'ipip)sub calculated with HISQ/tree action and 
compared to the renormalized Polyakov loop and light quark number fluctuation 
discussed in relation to the deconfining transition. Interestingly, the rapid decrease 
in the subtracted chiral condensate happens at temperatures that are smaller than the 
temperatures where the Polyakov loop rises rapidly as first noticed in Ref . [52] . On the 
other hand the rapid change in xi and {ip'ip)sub happen roughly in the same temperature 
region. The above results for {ipip)sub obtained with HISQ/tree action agree well with 
the continuum extrapolated stout results. 
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Figure 11. The subtracted chiral condensate calculated with HISQ/tree action on 
Nr = 8 lattice and compared to the renormalized Polyakov loop (left) and light quark 
number fluctuation (right). 



Another way to get rid of the niultiphcative and additive renormahzation is to 
subtract the zero temperature condensate and multiply the difference by the strange 
quark mass, i.e. consider the following quantity 



A 



R 



d + nfmsr\{{ipi))q^^ - {i'i')g,o), 



q = l,s. (36) 

1 for the strange quark, while d 



As before, nj = 2 for the light quarks and Uf 
is a normalization constant. The factor rf was introduced to make the combination 
dimensionless. It is convenient to choose the normalization constant to be the light 
quark condensate for mi = multiplied by rrigrf. In Fig. [12] the renormalized quark 
condensate is shown as function of the temperature for HISQ/tree and stout actions. 
We see a crossover behavior for temperature of (150 — 160)MeV, where Af drops by 
50%. The difference between the stout and HISQ/tree results is a quark mass effect. 
Calculations for HISQ/tree action were performed for mi = ms/20, while the stout 
calculations were done for the physical light quark masses, mi = ms/27.3. For a direct 
comparison with stout results, we need to extrapolate the HISQ/tree data in the light 
quark mass and also take care of the residual cutoff dependence in the HISQ/tree data. 
This was done in Ref. [56| and the results are shown in the figure as black diamonds 
demonstrating a good agreement between HISQ/tree and stout results. Contrary to 
Af the renormalized strange quark condensate Af shows only a gradual decrease over 
a wide temperature interval dropping by 50% only at significantly higher temperatures 
of about 190 MeV. 



5.2. The chiral susceptibility 

For a true chiral phase transition the chiral susceptibility diverges at the transition 
temperature. For physical value of the quark masses we expect to see a peak in the 
chiral susceptibility at a certain temperature that defines the crossover temperature. The 
chiral susceptibility also needs a multiplicative and additive renormahzation. Therefore 
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Figure 12. The renormalized chiral condensate A^ for the HISQ/tree actfon 
and mi/nis — 0.05 compared to the stout data. In the right panel, we show the 
renormahzed strange quark condensate Af for the HISQ/tree action. 



45 
40 
35 
30 
25 
20 
15 
10 
5 




Xr/T^ 




HISQ/tree: N.,=6 * 

N.,=8m 

asqtad: N^=8 @ 

N,=12« 

*^ stout: N.j=8 A 

N=10 A 



% h ig [MeV] 



120 140 160 180 200 220 240 




HISQ/tree: N^=6 * 

N,=8m 

asqtad: N^=8 ® 

N.,=12« 

stout: N^=8 A 

N =10a 



, T [MeV] 
220 240 



Figure 13. The renormahzed chiral susceptibility XR for the asqtad and HISQ/tree 
actions obtained at mi — 0.05tos and compared with the stout action results [52]. The 
temperature scale for HISQ/tree and asqtad actions is set using ri (/k) in the left 
(right) panels. 



the following quantity is considered 



(Xm,/(^) -XmAT = 0)). 



(37) 



2^4 2^4 

The numerical results for this quantity are shown in Fig. [13] for HISQ/tree, asqtad 
and stout actions Q There is a fairly good agreement between the results obtained with 
different actions if fx scale is used. At low temperatures the stout results are above the 
asqtad and HISQ/tree results due to quark mass effects [56] . 

* In Ref. [F5 the light quark mass was used instead of m^ in Eq. (|57]). In the comparison this was 
taken into account [55]. 
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5.3. 0(N) scaling and the transition temperature 

As discussed in the beginning of this section for physical nig and vanishing hght quark 
masses the transition is expected to be second order belonging to the 0(4) universality 
class in the continuum limit. At non-zero lattices spacing, however, we should expect 
0(2) universality class as only a part of the chiral symmetry is preserved in the staggered 
fermion formulation. Therefore in the following we will use the term 0{N) universality 
class as referring to either 0(4) or 0(2) universality class. In the vicinity of the chiral 
phase transition, the free energy density may be expressed as a sum of a singular and a 
regular part, 

T 

f = -—InZ = fsingit, h) + freg{T, Till, 171,) . (38) 

Here t and h are dimensionless couplings that control deviations from criticality. They 
are related to the temperature T and the light quark mass mi as 

i = ilzi? , ft= Iff , H^^. (39) 

to T^ ho nis 

where T^ denotes the chiral phase transition temperature, i.e., the transition 
temperature at iJ = 0. The scaling variables t, h are normalized by two parameters to 
and /lo, which are unique to QCD and similar to the low energy constants in the chiral 
Lagrangian. These need to be determined together with T°. In the continuum limit, 
all three parameters are uniquely defined, but depend on the value of the strange quark 
mass. On the lattice they will depend on N^ as well. 

The singular contribution to the free energy density is a homogeneous function of 
the two variables t and h. Its invariance under scale transformations can be used to 
express it in terms of a single scaling variable 



to TO \H) zo TO \H. 

where (5 and 5 are the critical exponents of the 0{N) universality class and zq = to/Z^o • 
Thus, the dimensionless free energy density / = //T*^ can be written as 

f{T,mi,m,) = h^^^"ff{z) + freg{T,H,ms) , (41) 

where // is the universal scaling function and the regular term f^eg gives rise to scaling 
violations. This regular term can be expanded in a Taylor series around (t, h) = (0, 0). 
It should be noted that the reduced temperature t may depend on other parameters 
which do not explicitly break chiral symmetry. In particular, it depends on light and 
strange quark chemical potentials Hg, which in leading order enter only quadratically, 

^"'■° + i:n(^)V-..^^). (42) 

q=l,s ^ ^ / 
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The transition temperature can be defined as peaks in the susceptibihties (response 
functions) that are second derivatives of the free energy density with respect to relevant 
parameters. Since there are two relevant parameters we can define three susceptibilities: 

Thus, three different pseudo-critical temperatures T^,/, Tti and T^j can be defined. 
In the vicinity of the critical point the behavior of these susceptibilities is controlled 
by three universal scaling functions that can be derived from ff. In the chiral limit 
Tm^i = Tt^i = Tt^t = T^- There is, however, an additional complication for 0{N) 
universality class: while Xm,i and Xt,i diverge at the critical point for m; — )• 

Xm,i^rni' , Xt,i^m\ , (44) 

Xt,t is finite because a < for 0{N) models ( Xt,t ~ \A~" )■ Therefore, one has to 
consider the third derivative of / with respect to t : 

Xtu = -g^- (45) 

In the vicinity of the critical point the derivatives with respect to t can be estimated 
by taking the derivatives with respect to ^f, i.e. the response functions Xt,i and Xt,t,t 
are identical to the second Taylor expansion coefficient of the quark condensate and the 
sixth order expansion coefficient to the pressure, respectively. The former controls the 
curvature of the transition temperature as function of the quark chemical potential fiq 
and was studied for p4 action using A',- = 4 and 8 lattices |142j . The later corresponds to 
the sixth order quark number fluctuation which is related to the deconfinement aspects 
of the transition. The fact that this quantity is sensitive to the chiral dynamics points 
to a relation between deconfining and chiral aspects of the transition. In the following 
I discuss the determination of the transition temperature defined as peak position of 

Xm,h l-^- -^ c -'■m,l- 

5.4- Determination of the transition temperature 

The 0{N) scaling described in the above subsection can be used to determine the 
pseudo-critical temperature of the chiral transition. For the study of the 0{N) scaling 
it is convenient to consider the dimensionless order parameter 

Mb = rris^^. (46) 

The subscript "b" refers to the fact that this is a bare quantity since the additive UV 
divergence is not removed. From the point of view of the scaling analysis this divergent 
term is just a regular contribution. For sufficiently small quark mass and in the vicinity 
of the transition region we can write 

M,{T, H) = h^/'fcit/h"^') + fM,reg{T, H) . (47) 
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Here faiz) is the scaling function related to // and was calculated for 0(2) and 0(4) 
spin models |143[ I144[ I145[ I146[ 1147] . The regular contribution can be parametrized as 
[56] 



fM,reg{T,H)=at{T)H 



Oo + ctr 



m TiO 



+ a2 



T 




C \ C 

Then we have the following behavior for the light chiral susceptibility 

Til 

hn 



(4^ 



Xm,l 



J^2 



mt 



-1/^1/5-1 f (^\ , dfM,reg{T,H) \ 
h ^^^ ^ OH J 



with f^{z) 



][fG{z) - ^fciz)]. 



(49) 



One performs a simultaneous fit to the lattice data for Mf, and Xm,i treating 
T°, to, ^05 O'O, cti and 02 as fit parameters [56]. This gives a good description of the quark 
mass and temperature dependence of Xm,i and allows to determine accurately the peak 
position in Xm,i- As an example in Fig. I show the 0(4) scaling fits for iV^ = 8 lattice 
data obtained with HI SQ /tree action. The scaling fit works quite well. Similar results 
have been obtained for A^.^ = 6 and 12 as well as for asqtad action on A^,- = 8 and 12 
lattices [56]. Furthermore, scaling fits have been performed assuming 0(2) universality 
class. The quality of these fits were similar to the 0(4) ones and the resulting transition 
temperatures turned out to be the same within statistical errors [56] . Having determined 
Tc for HISQ/tree and asqtad action for each Nr a combined continuum extrapolation 
was performed using different assumption about the A^,- dependence of T^. This analysis 
resulted in [56l: 



Te = (154 ± 9)MeV. 
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Figure 14. Scaling fits and data for the chiral condensate Mb calculated with 
the HISQ/tree action on lattices with temporal extent Nr = 8 (left) and the chiral 
susceptibility Xm,i(i'ight). The data for Mb at mi/nis — 0.025 and for Mb and Xmj at 
mi/ms = 0.05 are fit simultaneously using the 0(4) scaling Ansatz. The points used in 
the scaling fits are plotted using open symbols. The dotted lines give the data scaled 
to the physical quark masses. 
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The analysis also demonstrated that HISQ/tree and asqtad action give consistent 
results in the continuum limit. The Budapest- Wuppertal collaboration found T^ = 
147(2) (3)MeV, 157(3) (3)MeV and 155(3) (3)MeV defined as peak position in xr, 
infiection points in {4'4')sub and Af respectively [71]. These agree with the above value 
within errors. 

5. 5. The chiral transition temperature for small chemical potential 

The scaling analysis described in the above subsection is also useful for the determination 
of the chiral transition temperature for small quark chemical potential. Also for non- 
zero quark chemical potential the phase transition happens at t = 0. Together with Eq. 
P2|) this implies that the fig dependence of the chiral phase transition temperature is 
given by 

TeV) = T»-«:.feV + 0|(^) 1. (51) 




The value of k/ can be determined through lattice calculations of the mixed susceptibility 
Xt,i [142] . The mixed susceptibility can be expressed in terms of the derivative of the 
scaling function of the order parameter fciz) 

^ = ^-^h-^'-^y^^y^z). (52) 

T torris 
Calculations with p4 action on A^,- = 4 and 8 lattices of the mixed susceptibility have 
been performed for several quark masses and resulted in the value [142] 

Kz = 0.059(2)(4). (53) 

This result for the curvature of the critical line is about factor of two larger than the 
result obtained in Ref. |148[ 1149] using multi-parameter re-weighting technique. On 
the other hand it is consistent with the results obtained using the imaginary chemical 
potential technique |l50l [T5niT52] . 

Finally let me mention that the study of the mixed susceptibility allows to estimate 
the width of the chiral crossover for physical light quark masses. The pseudo-critical 
temperatures Tt^i and Tm,i are determined by the peak positions Zp of the fciz) and 
fx{z) when scaling violations are neglected. Namely, 

lElzll='jLH'm), p = m,t. (54) 

The peak positions of the above scaling functions for 0(4) model have been determined 
to be Zp = 0.74(4) and 1.374(30) respectively |147] . From these and from the values 
of zq and T° published in Ref. |142] for p4 action one can estimate Tm.i to be higher 
than Tt^i by 4MeV and 2MeV for physical value of the light quark mass for N^- = 8 
and 4 respectively. These estimates are reliable since the scaling violations were found 
to be small for p4 action |142] . In particular, the combination —mstoh^^~^y^^^^xt,i/T'^ 
as function of z has a peak at 2; ~ 0.74 |142] . It remains to be seen how this picture 
changes when HISQ/tree action is used. 
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6. Equation of State 

The equation of state has been calculated using different improved staggered fermion 
actions p4, asqtad, stout and HISQ/tree. The calculation of thermodynamic 
observables proceeds through the calculation of the trace of the energy momentum 
tensor, e — 3p, also known as trace anomaly or interaction measure. This is due to the 
fact that this quantity can be expressed in terms of expectation values of local gluonic 
and fermionic operators, (see e.g. Ref. |122] ). Different thermodynamic observables can 
be obtained from the interaction measure through integration of the trace anomaly IJ. 
The pressure can be written as 

The lower integration limit Tq is chosen such that the pressure is exponentially small 
there. Furthermore, the entropy density can be written as s = (e + p)/T. Since the 
interaction measure is the basic thermodynamic observable in the lattice calculations it 
is worth to discuss its properties more in detail. In Fig. [15] (left panel) I show the results 
of calculation with pA and asqtad actions using A^,- = 6 and 8 lattices and light quark 
masses mi = rris/lO, where rris is the physical strange quark mass. These calculations 
correspond in the continuum limit to the pion mass of 220MeV and 260MeV for p4 and 
asqtad respectively. The interaction measure shows a rapid rise in the transition region 
and after reaching a peak at temperatures of about 200MeV decreases. Cutoff effects 
(i.e. Nr dependence) appears to be the strongest around the peak region and decrease 
at high temperatures. For temperatures T < 270MeV calculations with p4 and asqtad 
actions have been extended to smaller quark masses, mi = ms/20, that correspond to 
the pion mass of about IGOMeV in the continuum limit |1541I155] . It turns out that the 
quark mass dependence is negligible for mi < rris/lO. Furthermore, for astqad action 
calculations have been extended to Nr = 12 lattices |155l 1156] . The trace anomaly 
was calculated with HISQ/tree action on lattices with temporal extent N^- = 6 and 8 
and mi = ms/20 |155] 1156] (corresponding to ?7i^ = 160MeV in the continuum limit). 
Finally, calculation of the trace anomaly and the equation of state was performed with 
stout action using N^. = 4, 6, 8, 10 and 12 and physical hght quark masses |153] . 
Using the lattice data from A^^ = 6, 8, 10 a continuum estimate for different quantities 
was given |153] . In Fig. [15] (right panel) the results of different lattice calculations of 
e — 3p corresponding to the pion masses close to the physical value are summarized. 
I also compare the lattice results with the parametrization s95p-vl of e — 3p. This 
parametrization combines lattice QCD results of Refs. [871 1122] at high temperatures 
with hadron resonance gas model (HRG) at low temperatures (T < 170MeV) [135] . 
At low temperatures there is a fair agreement between the results obtained with stout 
action and the results obtained with HISQ/tree action as well as with asqtad action for 
Nr = 12. All these lattice results are slightly above the HRG curve. For N^- = 8 cutoff 
effects are significant for asqtad and p4 action. As the result the corresponding lattice 

jj A somewhat different approach was used in Ref. [153j 
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data fall below the HRG (s95p-vl) curve at low temperatures. When cutoff effects 
in the hadron spectrum are taken into account in the HRG model a good agreement 
between the p4 data and HRG result can be achieved |135] . The height of the peak in 
e — 3p is the same for asqtad and HISQ/tree actions. At the same time it is smaller 
for the stout action. Since the dominant cutoff effects of order a'^T'^ are eliminated, the 
A^T--dependence is expected to be small for p4, asqtad and HISQ/tree action at high 
temperatures. We see that for T > 250MeV all these lattice actions lead to similar 
results. At temperatures above 350MeV we also see a good agreement with the stout 
results. 
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Figure 15. The interaction measure calculated with mi = ms/10 and p4 and asqtad 
actions [122) (left) and with nii = ms/20 ior pA action |154j as well as with HISQ/tree 
and asqtad actions |155| . Also shown in the figure are the continuum estimates 
obtained with stout action and the parametrization based on hadron resonance gas 
(HRG) model [T35] . 



The pressure, the energy density and the entropy density are shown in Fig. [161 The 
energy density shows a rapid rise in the temperature region (170 — 200) MeV and quickly 
approaches about 90% of the ideal gas value. The pressure rises less rapidly but at the 
highest temperature it is also only about 15% below the ideal gas value. In the previous 
calculations with the p4 action it was found that the pressure and the energy density are 
below the ideal gas value by about 25% at high temperatures |157j . A possible reason 
for this larger deviation could be the fact that the quark masses used in this calculation 
were fixed in units of temperature instead being tuned to give constant meson masses 
as lattice spacing is decreased. As discussed in Ref. |158] this could reduce the pressure 
by 10 — 15% at high temperatures. In Fig. [16] I also show the entropy density divided 
by the corresponding ideal gas value and compare the results of lattice calculations 
with resummed perturbative calculation |159[ 1160] as well as with the predictions from 
AdS/CFT correspondence for the strongly coupled supersymmetric Yang-Mills theory 
|161] . The later is considerably below the lattice results. Note that the pressure, the 
energy density and the trace anomaly have also been recently discussed in the framework 
of resummed perturbative calculations which seem to agree with lattice data quite well 
at high temperatures |162] . 
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The difference between the stout action and the pA and asqtad actions for the 
trace anomaly translates into the differences in the pressure and the energy density. 
In particular, the energy density is about 20% below the ideal gas limit for the stout 
action. 
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Figure 16. The energy density and tlie pressure as function of the temperature (left), 
and the entropy density divided by the corresponding ideal gas value (right). The 
dashed lines in the right panel correspond to the resummed perturbative calculations 
while the solid black line is the AdS/CFT result. 



7. Meson correlators and spectral functions 



In-medium meson properties as well as some transport coefficients are encoded in meson 
spectral functions. Medium modification of meson spectral functions can serve as 
diagnostic tools of the medium created in heavy ion collisions. For example, the excess 
in the low mass dilepton rate can be related to the modification of the light vector meson 
spectral function (see Ref. |163] for a recent review). The suppression of quarkonium 
yield was suggested by Matsui and Satz as signal for QGP formation in heavy ion 
collisions |164] . Studying quarkonium melting in QGP and the consequent suppression 
of the quarkonium yield is subject of a large experimental and theoretical effort (see 
Refs. [1651 [1661 MH [168] for reviews). 

The spectral function aH{u},p) for a given meson channel if in a system at 
temperature T can be defined through the Fourier transform of the real time two point 
functions D^ and D^ or equivalently as the imaginary part of the Fourier transformed 
retarded correlation function [57] . 
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Table 1. Meson states in different channels for light, charm and bottom quarks. 



In essence an is the Fourier transformation of the thermal average of the commutator 
[J(a:),J(0)]. 

In the present paper we study local meson operators of the form 

Jnit, x) = q{t, x)THq{t, x) (58) 

with q{t,x) is the quark field operator and 

r^ = 1, 75, 7a«> 757m' 7a.7i^ (59) 

for scalar, pseudo-scalar, vector, axial- vector and tensor channels. The relation of these 
quantum number channels to different meson states is given in Tab. [H In the vector 
channel I will use the subscript V to denote the sum over all four components, while 
the subscript ii will be used for the spatial components of the vector correlators and 
spectral functions. 

The correlators D^^ {t,x) satisfy the well-known Kubo-Martin-Schwinger (KMS) 
condition |57] 

D>{t,x) = D<{t + i/T,x). (60) 

Inserting a complete set of states and using Eq. (160|) . one gets the expansion 

c^h(w,p) = 



^e-^"/^x 



{n\MO)\m)\' {S\p, + k; 



^D 



5\p, + k 



fc;; 



(61) 



where Z is the partition function, /c"(™) refers to the four-momenta of the state \n{m)) 
andp^ = {uj,p). 

A stable meson state contributes a S function-like peak to the spectral function: 

aniuj^p) = mJH\H)\hiu)6ip' - M^), (62) 

where Mh is the mass of the state and e(po) is the sign function. For a quasi-particle in 
the medium one gets a smeared peak, with the width being the thermal width. As one 
increases the temperature the width increases and at sufficiently high temperatures, the 
contribution from the meson state in the spectral function may be sufficiently broad. 
At some point it is not very meaningful to speak of it as a well defined state any more. 
The spectral function as defined in Eq. (16T1) can be directly accessible by high energy 
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heavy ion experiments. More precisely, the spectral function for the vector current is 
directly related to the differential thermal cross section for the production of dilepton 
pairs [THQtITTn]: 



dW 



p=0 ^ ' 



duj(F'p 

Then presence or absence of a bound state in the spectral function will manifest itself 
in the peak structure of the differential dilepton rate. 

In finite temperature lattice calculations, one calculates Euclidean time 
propagators, usually projected to a given spatial momentum: 

Guir^p) = Jd^xe'P^-^{T,JH{T,x)MO,0)) (64) 

This quantity is an analytical continuation of D-^[xq,p) 

GH{r,p) = D>{-zT,p). (65) 

The KMS condition implies the following relation for the Fourier transforms of D^ and 

D<{uj,p) = D>i-co,p) = e--/'^D>{co,p). (66) 

This leads to a relation between D^{uj,p) and the spectral function 

D>{co,p) = 27raH{oo,p)e-/^/{e^/^ - 1). (67) 

Using this and Eq. (!65|) we get the following integral representation for the Euclidean 
time correlator 



roc 
Gh{t,p)= / dujGH{uJ,p)K{(XI,T) 

Jo 



cosh(a;(r-l/2T)) 
K[uj,T) = — -— — . (o8) 

This equation is the basic equation for extracting the spectral function from meson 
correlators. Equation (1681) is valid in the continuum. Formally the same spectral 
representation can be written for the Euclidean correlator calculated on the lattice 
G^*(r, jj). The corresponding spectral function, however, will be distorted by the effect 
of the finite lattice spacing, in particular, the spectral function is zero above certain 
energy u > Umax- These distortions have been calculated in the free theory |171l I172J . 
When discussing the numerical results in following sections the subscript H denoting 
different channels for meson correlators and spectral functions will be omitted unless 
stated otherwise. 

To get some information on the spectral functions from lattice QCD the 
corresponding Euclidean time correlation functions have to be calculated at as many 
separations in the time direction as possible, i.e. one should use large N^-. One way 
to accomplish this task with available computer resources is to use anisotropic lattices, 
i.e. lattices with different spacings at and a^ in time and space directions, such that 
^ = at/ as > 1. In addition one also often uses the quenched approximation for studying 
spectral functions on the lattice. 
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The low frequency limit of the vector spectral function gives information about the 
transport coefficients of the medium, namely the electric conductivity and the heavy 
quark diffusion constant. Other transport coefficients, such as the shear and bulk 
viscosity are related to the correlation functions of gluonic operators (see Ref. |173] 
for a review). The quark flavor diffusion constant Vq can be deflned through the time 
derivative of the quark number density in the rest frame of the thermal system as follows 

dtnq = VqV\ + 0{V'). (69) 

This equation holds for small deviation from thermal equilibrium. The quark flavor 
diffusion constant can be determined from the spatial component of the vector spectral 
function as follows 

V,^J-Mn.^-^, (70) 

6Xq ^-^0 UJ 

here Xq is the quark number susceptibility. The electric conductivity is related to the 
correlation function of electric currents 

Jr = j:QMi^^hMt,^). (71) 

/ 

with Qf being the electric charge of quark flavor /. The correlator of the electromagnetic 
current receives contribution from quark line connected and quark line disconnected 
diagrams. At high temperatures the contribution from the disconnected diagrams is 
small. For degenerate quark flavors the electric conductivity can be written as (174j 



^q ^q,con ~r LJq^dis' v'-^J 



For degenerate m, d and s quarks the disconnected contribution vanishes since 
Y.f=u,d,s Qf = 0- The quark flavor diffusion coefficient for single heavy quark flavor 
{ruq ^ T) is called the heavy quark diffusion constant D. It is related to the momentum 
drag coefficient tj in the Langevin dynamics of heavy quarks |175] : 

T 

D = -—. (73) 

mgT] 

In lattice QCD one also calculates meson correlation function in one of the spatial 
directions, say z 

G{z, T)= I dxdy f dT{J{x, y, z, r) J(0, 0, 0, 0)). (74) 

The spatial correlation function is related to the meson spectral function at non-zero 
spatial momentum 

G{z,T) = r dp^e'P^^ /■°°^^ cT(a;,p„T) ^ ^^^^ 

J-oo Jo UJ 

Thus the temperature dependence of the spatial correlation function also provides 
information about the temperature dependence of the spectral function. Medium effects 
are expected to be the largest at distances which are larger than 1/T. At these distances 
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G{z, T) decays exponentially and this exponential decay is governed by a screening mass 
Mgcr- If there is a lowest lying meson state of mass M, i.e. the spectral function can 
be well approximated by Eq. f l62l) . then the long distance behavior of the spatial meson 
correlation function is determined by the meson mass, i.e. Mgcr = M. At very high 
temperatures the quark and anti-quark are not bound and the meson screening mass is 



given by 2w (TrT)^ + m^, where m^ is the quark mass and ttT is the lowest Matsubara 
frequency. Therefore, a detailed study of spatial meson correlators and screening masses 
can provide some information about the melting of meson states at high temperatures. 
One would like to obtain the spectral functions through lattice calculation of 
the temporal meson correlation function. The obvious difficulty in the reconstruction 
of the spectral function from Eq. f l68|) is the fact that the Euclidean correlator is 
calculated only at (9(10) data points on the lattice, while for a reasonable discretization 
of the integral in Eq. ( I68l) we need (9(100) degrees of freedom. The problem can be 
solved using Bayesian analysis of the correlator, where one looks for a spectral function 
which maximizes the conditional probability P[a\DH] of having the spectral function 
a given the data D and some prior knowledge H (for reviews see |176[ 1177] ). Different 
Bayesian methods differ in the choice of the prior knowledge. One version of this 
analysis which is extensively used in the literature is the Maximum Entropy Method 
(MEM) |178[ 1179] . It has been used to study different correlation functions in QCD 
[1791 [T76l[l80l[T8ll[T82l[l83l[l8l[T85lE86l[l8a [193]. In this 

method the basic prior knowledge is the positivity of the spectral function and the prior 
knowledge is given by the Shannon - Janes entropy 



S = i duo 



oiuo) — niiuo) — (yiuj) ln( — — — ) 

m[uj) 



(76) 



The real function m{uj) is called the default model and parametrizes all additional prior 
knowledge about the spectral functions, e.g. such as the asymptotic behavior at high 
energy |179[ 1176] . For MEM the conditional probability can be written as 

P[a\DH] = eM-\x' + o^S), (77) 

with x^ being the standard likelihood function and a is a real parameter. 

Early results on meson spectral functions obtained with MEM gave a number of 
unexpected results. In particular, the MEM analysis of the charmonium correlators 
calculated in quenched QCD suggested that IS charmonium states can survive in QGP 
up to temperatures as high as 1.6Tc (see e.g. )188j ). Here T^ is the deconfinement 
phase transition temperature of SU{'i) gauge theory. This is in odds with the expected 
medium modification of the heavy quark potential that in addition to the effects of 
color screening also includes an imaginary part [73| 1194] . In fact, many potential 
model calculations indicate melting of IS charmonium state in the deconfined medium 
|195i I196i I197i I198[ I199i 1200] . Furthermore, even in the light meson spectral functions 
peak structures have been observed [185] . It is not clear to what extent these correspond 
to physical effects or are artifacts of the calculations. Next we will consider the 
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Figure 17. The spectral function in the pseudo-scalar channel calculated for 
several lattice spacings (left). The default model dependence of the spectral function 
calculated on 24^ x 160 lattice (right). The inset in the right panel shows the different 
default models used in the analysis. 



determination of the meson spectral function separately for the case of heavy quarks 
and light quarks. 



7. 1 . Quarkonium spectral functions at zero temperature 

The best way to understand the structure of meson spectral functions and the issues 
related to their determination on the lattice is to consider the case of heavy quarkonium 
at zero temperature. Charmonium correlation functions have been studied in detail in 
Ref. |201j using quenched anisotropic lattices and the so-called Fermilab formulations 
for heavy quarks |202] . Spectral functions have been extracted using MEM. The spectral 
functions in the pseudo-scalar channel calculated at three lattice spacings are shown in 
Fig. [TTl The first peak corresponds to the ground state, ric{lS). The second peak is 
actually the combination of several excited states, since MEM cannot resolve individual 
excited meson states due to the small splitting between them |201] . For a; > 5 GeV we 
see the continuum. The spectral function becomes zero above certain energy u due to the 
presence of finite lattice as expected in the free theory. In Fig. [T7|I also show the spectral 
function calculated for different default models m{u). The default model dependence 
is small for u < 5GeV, while above that energy it is significant due to the fact that 
there are only very few data points that carry information about the spectral function 
in that region. As the temperature increases the number of data points available for 
the analysis as well as the maximal extent of the time direction r^ax = 1/(2T) become 
smaller. To demonstrate this point in Fig. [18] I show the zero temperature charmonium 
spectral function calculated with different number of data points and T^ax that are 
characteristic for temperatures corresponding to QGP. As Tmax decreases MEM looses 
the ability to reconstruct the peak corresponding the excited states and the ground state 
peak is significantly broadens. For Tmax = 0.22fm even the position of the ground state 
peak is not reproduced correctly. Furthermore, the shape of the spectral functions, in 
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Figure 18. The dependence of the pseudo-scalar spectral function on t„ 



particular, the presence of the ground state peak also becomes sensitive to the default 
model. 

Since the analysis of the quarkonium spectral functions using MEM turns out to be 
quite complicated already in the zero temperature limit, it is important to understand 
the temperature dependence of the quarkonium correlation functions and try to identify 
possible sources of the temperature dependence that are related to melting of the bound 
states. From Eq. (!68|) it is clear that the temperature dependence of the meson 
correlation functions comes from two sources: the trivial temperature dependence of 
the integration kernel K{t, u) and the temperature dependence of the spectral function 
a{u,T). To get rid of the first trivial temperature dependence one can consider the 
reconstructed correlation function 



("OO 

GreciT,T) = / dua{u,T = 0)K{u,t). 
Jo 



(7J 



If the spectral function does not change across the deconfinement transition 
G{T,T)/Grec{'^,T) should be unity. Deviations of this ratio from unity indicate 
temperature dependence of the spectral functions. The ratio G{T,T)/Grec{'T,T) was 
first studied in Ref. [189] and subsequently in Refs. [HQl |193l [2011 1203]. It was 
found that in the pseudo-scalar channel this ratio stays close to one and shows only 
small temperature dependence. This seemed to support the conclusion based on the 
spectral functions extracted from MEM that ground state charmonium survives in the 
deconfined phase up to temperatures l.GTc. In the scalar and axial- vector channels large 
temperature dependence in G{T,T)/Grec{'T,T) was observed. This was interpreted as 
melting of IP charmonium states and fitted well into the expected picture of sequential 
melting. However, this conclusion was premature. In the deconfined phase there is an 
additional contribution to the spectral functions at very low frequency. This can be 
easily seen by calculating the spectral function in the free theory. In addition to the 
contribution to the spectral function that starts at twice the heavy quark mass there 
is a contribution proportional to uj6{oj) in all but the pseudo-scalar channel [172] . The 
delta function is smeared once the interactions of the heavy quarks with the medium 
are taken into account. In particular, in the vector channel, where this low u structure 
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Figure 19. The schematic structure of the vector spectral function for heavy quarks. 
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Figure 20. The ratio of the derivatives of the charmonium correlators to the 
corresponding reconstructed correlators for the pseudo-scalar (left) and scalar (right) 
channels |205j . 



is related to the heavy quark diffusion we have [204] 

1 rju 



u5{u) 



rj = T /[Dm, 



q>i 



(79) 



-K uj"^ + rf' 

with mq being the heavy quark mass and D being the heavy quark diffusion constant. 
The schematic structure of the spectral functions is shown in Fig. [191 The area under 
the transport peak is given by the quark number susceptibility Xq- For large quark 
mass the transport peak is narrow and is well separated from the high energy part of 
the spectral function that corresponds to bound states and/or unbound heavy quark 
anti-quark pairs. We expect similar structure in the spectral function in other channels 
as well. Thus we can write 

a{u, T) = (Tiow{uj, T) + ahigh{uj, T) (80) 

G{t,T) = Gion,{r,T) + GMgh{r,T). (81) 

Since the peak at low u is narrow we expect that the derivative of Giow{t,T) with 
respect of r is small, G'i^^{t,T) ^ 0. In Fig. [20] I show the ratio of the derivatives 
G'{t,T)/G'^^^{t,T) in the pseudo-scalar and scalar channels in the deconfined phase. 
As once can see from the figure these ratios show no strong temperature dependence 



Lattice QCD at non-zero temperature 



37 




10 15 

co[GeV] 




10 15 

ro[GeV] 



Figure 21. The spectral functions in pseudo-scalar channel calculated at T = 1.5Tc 
and compared to the zero temperature spectral function obtained using Tmax — 1/(2^)- 
In the left panel the spectral functions obtained with m(uj) = const are shown, while 
in the right panel the default model was obtained from the T — spectral function. 



and are close to one for all temperatures, including the highest temperature of 3Tc. 
This means that the large temperature dependence seen in G/Grec for the scalar and 
axial-vector channels comes mostly from the temperature dependence of Giow{t,T) 
and is not related to melting of P-wave charmonium states. A question arises how 
the observed temperature (in) dependence of Ghighij^T) is related to the expected in- 
medium modification of the heavy quark potential and melting of quarkonium states at 
sufficiently high temperatures. This question was addressed in Refs. |197] . where the 
quarkonium spectral functions have been calculated using potential model with screened 
potential. It was found that despite the significant change in the spectral functions the 
Euclidean correlator obtained from them do not change significantly in the deconfined 
medium and are compatible with the lattice results. Including the imaginary part in 
the analysis did not change this conclusion |199] . Attempts to reconstruct quarkonium 
spectral functions using MEM have been presented in Refs. [Ml UMl HHHl [HQl [20T] 
and no evidence for melting of quarkonium states has been found. However, in view of 
the difficulties related to the MEM analysis for small Tmax and limited number of data 
points one should be careful drawing conclusions. In particular the dependence of the 
results on the default model should be examined. In Fig. |21]I show charmonium spectral 
functions in the pseudo-scalar channel calculated at 1.5Tc and at zero temperature for 
two default models: m{ijj) = const and m{uj) that equals to the zero temperature 
spectral function for u > 5GeV and smoothly matched to a constant below that energy 
|201j . In both cases the difference between the spectral function at 1.5Tc and T = is 
small. However, the shape of the spectral function depends on the default model. While 
for m{u) = const we see a peak structure for more realistic choice of the default model 
this peak structure is absent even at zero temperature. Recent analysis of charmonium 
spectral function based on isotropic lattices that uses free lattice spectral function also 
does not find peak structure that can be associated with IS state in the deconfined phase 
|206] It is not clear to what extent the absence of bound state peaks is consequence of 
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quarkoniuin melting in QGP or due to the limited data set. It is clear, however, that 
existing MEM calculations of the spectral function do not provide evidence for existence 
of quarkonium bound states in QGP. 

The analysis described so far was performed in the quenched approximation. 
Calculation of the charmonium correlators and spectral functions was also performed 
in two flavor QCD |207j . The findings of this analysis are similar to the ones described 
above. Furthermore, in the case of the bottomonium correlators and spectral functions 
have been studied using NRQCD [SDHl EOH]. The main advantage of the NRQCD 
approach is that it allows to study correlators at larger Euclidean time separations, 
namely Tmax = 1/T. Furthermore, there is no low energy contribution to the spectral 
functions and thus the temperature dependence of the corresponding Euclidean time 
correlation function is directly related to the properties and melting of the bound 
states. A significant temperature dependence has been seen in the scalar bottomonium 
correlators that may indicate the melting of the P-wave bottomonium at temperatures 
slightly above the transition temperature. 

7.2. Meson spectral functions at non-zero temperature in the light quark sector 

Early attempts to calculate light meson spectral functions in QGP using MEM were 
presented in Refs. |180[ I183[ 11851 1187] . The results were inconclusive and to some 
extent confusing as peak structures in the spectral functions have been found up to 3Tc 
and no continuum was observed at high energies. It was pointed out, however, that even 
without using MEM Euclidean correlation functions can put stringent constraints on the 
spectral functions |180[ 1187] . In the vector channel the correlation function calculated 
on Nr = 12 and A^,- = 16 lattices was found to deviate from the free value by less than 
10% |180] . In the pseudo-scalar channel on the other hand a large enhancement over 
the free theory result was observed |187j possibly indicating non-perturbative effects in 
QGR 

More recently a detailed calculation of the vector correlation function on large 
quenched lattices was reported for T ~ 1.45Tc [210] . Calculations were performed using 
several lattice spacings. This enabled a reliable extrapolation to the continuum limit 
and it was found that the correlator never exceeds the free theory value by more than 
9% [210J. The corresponding results are shown in Fig. [22j Furthermore, the spectral 
functions reconstructed with MEM showed no evidence for peak like structures [210] . 
Since vector correlator is so close to the free limit the authors of Ref. [210 ] obtained the 
spectral function using the following model 

aii{u) = CBwXq o , /-n/o\2 + 7^(1 + ^)^'^ tanh(u;/4T)e(wo, A^), 

e{uo, A^) = (1 + e(-o-)/-^-)-i, (82) 

and treating cbw, T and k as fit parameters. Furthermore, several choices for the 
parameters uo and A^ have been considered, including uq = A^ = 0. These fits are 
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Figure 22. The vector correlation function, Gv{t,T) calculated on different lattices 
for T = 1.45rc (left) and the vector spectral function (right) obtained using the fit 
form given by Eq. (J82|) . 



shown in Fig. [82] and gave the following constrains for the electric conductivity [210] 

1 c 



1/3 < 



C T 

^--^ p.m. -*■ 



<1, Ce 



f 



(83) 



The disconnected part was neglected in this calculation. One can see from Fig. 
that while somewhat broad the transport peak is clearly visible in the vector spectral 
function contrary to the expectation based in strongly coupled supersymmetric Yang- 
Mills theory where the transport peak is absent |211j . On the other hand the analysis 
of the same lattice data using a different approach lead to much smaller value of the 
electric conductivity and no clear transport peak |174] . This may indicate that much 
more work is needed till a reliable result for the electric conductivity can be quoted. 



7. 3. Spatial meson correlation functions 

As discussed above spatial meson correlators are also sensitive to in-medium 
modification of the meson spectral functions. In particular, the change in the meson 
spectral functions is reflected in the meson screening masses, which should be close to 
the vacuum masses at low temperatures. At high temperatures, on the other hand, the 
meson screening masses should approach 27rT corresponding to the free quark limit. In 
the past meson screening masses have been studied in quenched approximation |212[ 1213] 
and the sharp change in their behavior was observed at the transition temperature. More 
recently spatial meson correlators and screening masses have been studied in 2-1-1 flavor 
QCD for physical strange quark mass and light quark masses mi = m,s/10 with p4 action 
|141] . In the case of full QCD rapid change of the spatial correlators and screening masses 
is expected in the vicinity of the chiral crossover. Numerical calculations show that this 
is indeed the case |141] . In Fig. [23] the meson screening masses in the pseudo-scalar 
and vector channels are shown as function of the temperature and the rapid change in 
the behavior of the screening masses takes place around 200MeV, which is close to the 
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Figure 23. Meson screening masses for light quarks in the pseudo-scalar and vector 
channels as function of the temperature calculated in 2+1 flavor QCD withp4 action on 
Nr = 6 (filled symbols) and Nr = 8 (open symbols) lattices |141| . The solid black line 
corresponds to the limit of free quarks 27rT. The dashed lines are the NLO results |214] 
evaluated with a 2-loop running coupling constant in Af 5-scheme for fi — ttT, 2t:T 
and AttT. 



chiral crossover temperature T^ for those lattices. The effective restoration of the chiral 
symmetry above the crossover temperature should manifest in approximate degeneracy 
of the vector and axial-vector correlation functions. Such degeneracy is observed at 
temperature 200MeV (c.f. Fig. 5 of Ref. |141j ). As discussed in section [5] the Ua{1) 
symmetry is expected to be restored at sufficiently high temperatures. The effective 
restoration of the Ua{^) should manifest itself in the approximate degeneracy of flavor 
non-singlet pseudo-scalar and scalar correlators. The difference between the scalar and 
pseudo-scalar screening masses rapidly decreases above the transition temperature but 
becomes compatible with zero only temperatures above 240MeV. Thus the effective 
restoration of the axial symmetry happens at T ~ 1.2Tc, with T^ being the chiral 
transition temperature. Calculations using DWF formulations confirm this result |140j . 

At high temperatures the screening masses can be calculated using weak coupling 
methods. The calculations of the screening masses up to next-to-leading order was 
performed in Ref. [214] and the first perturbative correction turned out to be positive, 
namely 

M,,, = 2nT + ^^g^il + Eo), (84) 

with Efj ~ 0.46939139 for Nf = 3 [214] . Furthermore, all the screening masses are 
expected to be degenerate at next-to- leading order |214j . In Fig. [23] we compare the 
lattice data in 2+1 fiavor case with the perturbative result of Ref. |214] . At high 
temperatures we still see significant differences between the pseudo-scalar and vector 
screening masses. Moreover the pseudo-scalar screening masses are smaller than 27cT. 
Note, however, that the lattice spacing (Nr) dependence of the screening masses is 
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non-negligible. Thus it will be important in the future to perform calculations on finer 
lattices to get control over the continuum limit. 

Spatial meson correlators and screening masses have been calculated also for 
charmonium |215] . The changes in the correlators and screening masses around the 
transition temperature are much smaller than for the light mesons. This fits into the 
picture that ground state charmonium does not melt at the crossover temperature. 
However, for T > 300MeV we see large change also for the charmonium screening 
masses. The behavior of the screening masses and their dependence on the spatial 
boundary conditions is compatible with the picture of unbound quarks |215j . This 
corroborates the discussion in the previous subsection that there is no evidence for 
survival of charmonium states at temperatures of l.ST^. 

7.4- Heavy quark diffusion constant 

Since for sufficiently heavy quarks the transport peak is very narrow it is very difficult if 
not impossible to determine the heavy quark diffusion using the procedure that was used 
for the electric conductivity qijj. Even for the smallest possible heavy quark diffusion 
constant D ~ l/(27rT) the low energy part of the vector correlator Giow(t, T) shows 
very little r-dependence [204] . An alternative way to estimate the heavy quark diffusion 
constant is integrate out the heavy quark fields and relate the drag coefficient to the 
correlation function of chromo-electric fields at time t and time connected by Wilson 
lines |217l 1218] . In particular, at leading order in the expansion in the inverse of heavy 
quark mass we have 



K = j dtGEit) 



GE{t) =^ I d'xiO^Eit, x)OU0, 0)) 

0%{t, x) = (p^t, x)gE\t, x)0(t, x) - 9\t, x)gE\t, x)e{t, x). (85) 

Here and 6 are the static quark and anti-quark fields, and E^ is the chromo-electric 
field. It is easy to generalize the above expression to Euclidean time and define the 
corresponding correlator Ge{t) which after integrating the heavy quark fields becomes 

The above correlation function was also calculated in perturbation theory to next-to- 
leading order |218] . The key difference between this correlation function and meson 
correlation function is that the transport contribution does not appear as a narrow 
peak in the corresponding spectral function This is principle should make it easier to 
estimate the corresponding transport coefficient. Numerical calculations of Ge{t) have 
been performed in SU{3) gauge theory on lattices with temporal extent A'^^ = 12 — 24 
|219[ 122m 1221] . The lattice data are significantly larger than the next-to- leading order 

ft attempts along these lines have been presented in Ref. [216] 
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result. Comparing the numerical results with a model correlation function based on 
perturbative calculations of Ref. |218] and strong coupling calculations in A^ = 4 Super- 
Yang-Mills theory |217] it was concluded that the heavy quark diffusion constant should 
be in the range 

D = {0.5- 1.0) /T. (87) 

The above value is considerably smaller than the perturbative estimate |175j and lies in 
the range used in phenomenological models [168j . 

8. Conclusions 

In recent years significant progress has been made in studying chiral and deconfining 
aspects of the QCD transition in finite temperature QCD. For several quantities that are 
relevant for the discussion of these aspects continuum extrapolation has been performed. 
Calculations by two groups that use different staggered quark formulation give results for 
the chiral transition temperature that agree in the continuum limit. It is now established 
that for physical value of chiral transition temperature is Tc ~ ISOMeV. Calculations 
using Wilson fermion or DWF formulation have been performed for physical or nearly 
physical quark masses and seem to confirm the staggered fermion results. The universal 
behavior of the chiral transition in the limit of vanishing light quark masses plays an 
important role also for the physical values of the light quark masses and allows to define 
the chiral crossover temperature in a meaningful way. 

The interplay between chiral and deconfinement aspects of the transition appears 
to be more complicated than earlier lattice studies suggested. There is no transition 
temperature that can be associated with the deconfining aspects of the transition for 
physical values of the light quark masses. Furthermore, the behavior of the Polyakov 
loop suggests that color screening sets in at temperatures that are higher than the chiral 
transition temperature. The deconfinement aspects of the QCD transition have been 
also studied in terms of fluctuations of conserved charges. In principle these fluctuations 
are are also sensitive to universal aspects of the chiral transition. However, for the 
lowest order, i.e. quadratic fluctuations the contribution coming from the regular 
part of the free energy density dominates. This implies in particular, that inflection 
points of the quark number susceptibilities cannot be used to define the transition 
temperature. Higher order fluctuations of conserved charges are more sensitive to the 
singular contribution and therefore are more suitable to study the interplay between 
chiral and deconfining aspects of the QCD transition at finite temperature. The 
restoration of the Ua{1) symmetry at high temperature has been also discussed and 
it was found that it is effectively restored at temperatures T ~ 1.2Tc. 

In the high temperature region ( T > 300 MeV) lattice results have been compared 
with the results obtained in weak coupling approaches, which seem to capture the 
qualitative features of the lattice data. In some cases we see a good agreement at 
quantitative level. 
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Finally I discussed lattice results on meson correlation functions. These are useful 
to study the fate of quarkonium states in QGP as well as for determination of some 
transport coefficients. While a lot of progress has been made in studying meson 
correlation functions and extracting the corresponding spectral functions more work 
is needed to reach definitive conclusions. 
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